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Models for the seds of infrared galaxies are reviewed, with particular application to 
hyperluminous infrared galaxies. 
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1 Ingredients for model SEDs 



The key ingredients for models for the spectral energy distributions (seds) of 
infrared galaxies are (1) A model for interstellar dust grains. Detailed 
models have been provided by Mathis et al (1977), Draine and Lee (1984), 
Rowan-Robinson (1986, 1992), Desert et al (1990), Siebenmorgen and Krugel 
(1992), Dwek (1998). (2) An assumed density distribution for the dust. 
Detailed physical density profiles have been developed by Yorke (1977) for HII 
regions, by Efstathiou and Rowan- Robinson (1991) for protostars and by Ef- 
tathiou et al (2000) for starbursts. However most authors have followed the 
approach of Rowan- Robinson (1980), who explored models with N(r) oc r~@. 
Silva and Danese (1998) have investigated exponential density profiles for the 
cirrus component. (3) The geometry of the dust. The most common as- 
sumption is spherical symmetry. Axi-symmetric models have been investigated 
by Efstathiou and Rowan- Robinson (1990, 1991, 1995), Pier and Krolik (1992), 
Granato et al (1994, 1997), Silva et al (1998). Rowan-Robinson (1995) gave an 
approach to modelling emission from an ensemble of discrete compact clouds. 
(4) A radiative transfer code The first accurate spherically symmetric ra- 
diative transfer code was by Rowan-Robinson (1980) and the first accurate 
axisymmetric code by Efstathiou et al (1990). Pier and Krolik (1992), Krugel 
and Siebenmorgen (1994), Granato et al (1997) and Silva et al (1998) have 
also developed axisymmetric codes. (5) Source components Typically, seds 
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can be understood as a superposition of several distinct source components. 
Rowan- Robinson and Crawford (1989) were able to fit the IRAS colours and 
spectral energy distributions of galaxies detected in all 4 IRAS bands with a 
mixture of 3 components, emission from interstellar dust ('cirrus'), a starburst 
and an AGN dust torus. Recently Xu et al (1998) have shown that the same 
3-component approach can be used to fit the ISO-SWS spectra of a large sam- 
ple of galaxies. To accomodate the Condon et al (1991) and Rowan-Robinson 
and Efstathiou (1993) evidence for higher optical depth starbursts, Ferrigno et 
al (2000, in prep.) have extended the Rowan-Robinson and Crawford (1989) 
analysis to include a fourth component, an Arp220-like, high optical depth 
starburst, for galaxies with log Lqq > 12. Efstathiou et al (2000) have given 
improved radiative transfer models for starbursts as a function of the age of 
the starburst, for a range of initial dust optical depths, and similar analysis 
has been carried out by Silva et al (1998) and Siebenmorgen (1999, 2000 - 
these proceedings). 

Sanders et al (1988) proposed, on the basis of spectroscopic arguments for a 
sample of 10 objects, that all ultraluminous infrared galaxies contain an AGN 
and that the far infrared emission is powered by this. Sanders et al (1989) pro- 
posed a specific model, in the context of a discussion of the infrared emission 
from PG quasars, that the far infrared emission comes from the outer parts of 
a warped disk surrounding the AGN. This is a difficult hypothesis to disprove, 
because if an arbitrary density distribution of dust is allowed at large distances 
from the AGN, then any far infrared spectral energy distribution could in fact 
be generated. Rowan-Robinson (2000a) discussed whether the AGN dust torus 
model of Rowan- Robinson (1995) can be extended naturally to explain the far 
infrared and submillimetre emission from hyperluminous infrared galaxies, but 
concluded that in many cases this does not give a satisfactory explanation. 
Rigopoulou et al (1996) observed a sample of ultraluminous infrared galaxies 
from the IRAS 5 Jy sample at submillimetre wavelengths, with the JCMT, 
and at X-ray wavelengths, with ROSAT. They found that most of the far in- 
frared and submillimetre spectra were fitted well with the starburst model of 
Rowan-Robinson and Efstathiou (1993). The ratio of bolometric luminosities 
at 1 keV and 60 /mi lie in the range 10~ 5 — 10~ 4 and are consistent with a 
starburst interpretation of the X-ray emission in almost all cases. Even more 
conclusively, Genzel et al (1998) have used ISO-SWS spectroscopy to show 
that the majority of ultraluminous ir galaxies are powered by a starburst 
rather than an AGN. 



2 Applications to quasars and hyperluminous infrared galaxies 

Quasars and Seyfert galaxies, on the other hand, tend to show a characteristic 
mid infrared continuum, broadly flat in vS v from 3 — 30/im. This component 
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was modelled by Rowan- Robinson and Crawford (1989) as dust in the narrow- 
line region of the AGN with a density distribution n(r) a r _1 . More realistic 
models of this component based on a toroidal geometry are given by Pier and 
Krolik (1992), Granato and Danese (1994), Efstathiou and Rowan-Robinson 
(1995). Rowan- Robinson (1995) suggested that most quasars contain both (far 
ir) starbursts and (mid ir) components due to (toroidal) dust in the narrow 
line region. 

One of the major discoveries of the IRAS mission was the existence of ultralu- 
minous infrared galaxies, galaxies with Lf ir > 10 12 h^QL Q (h 50 = H o /50). The 
peculiar Seyfert 2 galaxy Arp 220 was recognised as having an exceptional far 
infrared luminosity early in the mission (Soifer et al 1984). 

The conversion from far infrared luminosity to star formation rate has been 
discussed by many authors (eg Scoville and Young 1983, Rowan- Robinson et 
al 1997). Rowan- Robinson (2000b, in prep.) has given an updated estimate of 
how the star-formation rate can be derived from the far infrared luminosity, 
finding 

M^ all /[L m /L Q ] = 2.2 0/e xlO" 10 

where takes account of the uncertainty in the IMF (= 1, for a standard 
Salpeter function) and e is the fraction of uv light absorbed by dust, estimated 
to by 2/3 for starburst galaxies (Calzetti 1998). We see that the star-formation 
rates in ultraluminous galaxies are > 10 2 M Q yr~ 1 . However the time-scale of 
luminous starbursts may be in the range 10 7 — 10 8 yrs (Goldader et al 1997), 
so the total mass of stars formed in the episode may typically be only 10% of 
the mass of a galaxy. 

Here I discuss an even more extreme class of infrared galaxy, hyperluminous in- 
frared galaxies, which I define to be those with rest-frame infrared (1-1000 jum) 
luminosities, L bo i yir) in excess of lO 13 ' 22 /i5 O 2 L (=lO 13 '°/i6 5 2 L ). For a galaxy 
with an M82-like starburst spectrum this corresponds to Lqq > 1O 13 /i^ o 2 L , 
since the bolometric correction at 60 fim is 1.63. While the emission at rest- 
frame wavelengths 3-30 /xm in these galaxies is often due to an AGN dust torus 
(see below), I argue that their emission at rest-frame wavelengths > 50/xm is 
primarily due to extreme starbursts, implying star formation rates in excess of 
1000 M /yr. These then are excellent candidates for being primeval galaxies, 
galaxies undergoing a very major episode of star formation. More details are 
given in Rowan- Robinson (2000a). 

For a small number of these galaxies we have reasonably detailed continuum 
spectra from radio to uv wavelengths. Figures 1-6 show the infrared continua of 
some of these hyperluminous galaxies, with fits using radiative transfer models 
(specifically the standard M82-like starburst model and an Arp220-like high 
optical depth starburst model from Efstathiou et al (2000) and the standard 
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AGN dust torus model of Rowan- Robinson (1995). 
I now discuss the individual objects in turn: 
F10214 + 4724 

The continuum emission from F10214 was the subject of a detailed discussion 
by Rowan-Robinson et al (1993). Green and Rowan-Robinson (1995) have dis- 
cussed starburst and AGN dust tori models for F10214. Fig 1 shows M82-like 
and Arp 220-like starburst models fitted to the sumbillimetre data for this 
galaxy. The former gives a good fit to the latest data. The 60 jum flux requires 
an AGN dust torus component. To accomodate the upper limits at 10 and 
20 jum, it is necessary to modify the Rowan-Robinson (1995) AGN dust torus 
model so that the maximum temperature of the dust is 1000 K rather than 
1600 K. I have also shown the effect of allowing the dust torus to extend a 
further factor 3.5 in radius. This still does not account for the amplitude of 
the submm emission. The implied extent of the narrow-line region for this ex- 
tended AGN dust torus model, which we use for several other objects, would 
be 326 (L bol /10 13 L Q y/ 2 pc consistent with 60-600 (L^/IO 13 ^©) 1 / 2 pc quoted 
by Netzer (1990). Evidence for a strong starburst contribution to the ir emis- 
sion from F10214 is given by Kroker et al (1996) and is supported by the high 
gas mass detected via CO lines. Granato et al (1996) attempt to model the 
whole sed of F10214 with an AGN dust torus model, but still do not appear 
to be able to account for the 60 /zm emission. 

SMMJ02399 - 0136 

A starburst model fits the submm data well and the ISO detection at 15 
fim gives a very severe constraint on any AGN dust torus component. The 
starburst interpretation of the submm emission is supported by the gas mass 
estimated from CO detections (Frayer et al 1998). 

F08279 + 5255 

An M82-like starburst is a good fit to the submm data and an AGN dust torus 
model is a good fit to the 12-100 /xm data. The high gas mass detected via 
CO lines (Downes et al 1999) supports a starburst interpretation, though the 
ratio of L sb /M gas is on the high side. The submm data can also be modelled 
by an extension of the outer radius of the AGN dust torus. 

H1413 + 117 

The submm data is well fitted by an M82-like starburst and the gas mass im- 
plied by the CO detections (Barvainis et al 1994) supports this interpretation. 
The extended AGN dust torus model discussed above does not account for 
the submm emission. However Granato et al (1996) model the whole sed of 
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H1413 in terms of an AGN dust torus model. 
15307 + 325 

A starburst model gives a natural explanation for the 60-180 /zm excess com- 
pared to the AGN dust torus model required for the 6.7 and 14.3 /xm emission 
(Verma et al 2000), but the non-detection of CO poses a problem for a star- 
burst intepretation. 

PG1634 + 706 

The IRAS 12-100 fim data and the ISO 150-200 fim data (Haas et al 1998) 
are well-fitted by the extended AGN dust torus model and an upper limit can 
be placed on any starburst component. The non-detection of CO is consistent 
with this upper limit. 

In Fig 7 we show the far infrared luminosity derived for an assumed starburst 
component, versus redshift, for hyperluminous galaxies, with lines indicat- 
ing observational constraints at 60, 800 and 1250 /im. Three of the sources 
with (uncorrected) total bolometric luminosities above 10 14 /i^ 2 L o are strongly 
gravitationally lensed. IRAS F10214+4724 was found to be lensed with a mag- 
nification which ranges from 100 at optical wavelengths to 10 at far infrared 
wavelengths (Eisenhardt et al 1996, Green and Rowan- Robinson 1996). The 
'clover-leaf lensed system H1413+117 has been found to have a magnifica- 
tion of 10 (Yun et al 1997). Downes et al (1999) report a magnification of 
14 for F08279+5255. Also, Ivison et al (2000) estimate a magnification of 2.5 
for SMMJ02399-0136 and Frayer et al (1999) quote a magnification of 2.75 
±0.25 for SMMJ1401 1+0252. These magnifications have to be corrected for 
in estimating luminosities (and dust and gas masses) and these corrections are 
indicated in Fig 7. 

Although lensing may affect some of the most luminous objects, there is strong 
evidence for a population of galaxies with far ir luminosities (after correction) 
in the range l-3xl0 13 /i^ 2 L Q . I have argued that in most cases the rest-frame 
radiation longward of 50 /im comes from a starburst component. The luminosi- 
ties are such as to require star formation rates in the range 3-lOxlO 3 /i^ o 2 M yr, 
which would in turn generate most of the heavy elements in a 1O 11 M galaxy 
in 10 7 — 10 8 yrs. Most of these galaxies can therefore be considered to be 
undergoing their most significant episode of star formation, ie to be in the 
process of 'formation'. 

It appears to be significant that a large fraction of hyperluminous ir galaxies 
are Seyferts, radio-galaxies or QSOs. This is in part a selection effect, in that 
luminous AGN have been selected for submm photometric follow-up. For the 
12 objects found from direct follow-up of IRAS samples or 850 /mi surveys 
(Rowan- Robinson 2000a), 50 % are QSOs or Seyferts and 50 % are narrow- line 
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Fig. 1. Observed (rest-frame) spectral energy distribution for F10214, modelled with 
M82-type starburst (solid curve), Arp 220- type starburst (broken curve), AGN dust 
torus (dotted curve, modified AGN dust torus model - long-dashed curve). 




Fig. 2. Observed spectral energy distribution for SMMJ02399, notation as for Fig 
1. 

objects, similar to the proportions found for ultraluminous ir galaxies. How- 
ever despite the high proportion of ultraluminous and hyperluminous galaxies 
which contain AGN, this does not prove that an AGN is the source of the 
rest-frame far infrared radiation. The ISO-LWS mid-infrared spectroscopic 
programme of Genzel et al (1998) has shown that the far infrared radiation of 
most ultraluminous galaxies is powered by a starburst, despite the presence of 
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Fig. 3. Observed spectral energy distribution for F08279, notation as for Fig 1. 
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Fig. 4. Observed spectral energy distribution for H1413, notation as for Fig 1. 

an AGN in many cases. Wilman et al (1999) have shown that the X-ray emis- 
sion from several hyperluminous galaxies is too weak for them to be powered 
by a typical AGN. 

In the Sanders et al (1989) picture, the far infrared and submillimetre emission 
would simply come from the outer regions of a warped disk surrounding the 
AGN. Some weaknesses of this picture as an explanation of the far infrared 
emission from PG quasars have been highlighted by Rowan- Robinson (1995). 
A picture in which both a strong starburst and the AGN activity are triggered 
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F15307, z = 0.93 
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Fig. 5. Observed spectral energy distribution for F15307, notation as for Fig 1. 
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Fig. 6. Observed spectral energy distribution for PG1634, notation as for Fig 1. 

by the same interaction or merger event is far more likely to be capable of 
understanding all phenomena (cf Yamada 1994, Taniguchi et al 1999). 

Where hyperluminous galaxies are detected at rest-frame wavelengths in the 
range 3-30 /im (and this can correspond to observed wavelengths up to 150 
/im), the infrared spectrum is often found to correspond well to emission from 
a dust torus surrounding an AGN (eg Figs 1, 3, 4, 6). This emission often 
contributes a substantial fraction of the total infrared (1-1000 /im) bolometric 
luminosity. The advocacy of this paper for luminous starbursts relate only to 
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Fig. 7. Bolometric luminosity in starburst component for galaxies with luminosi- 
ties > 10 13 L Q (Tables 1-4). The galaxies labelled L are known to be lensed. Loci 
corresponding to the limits set by S(60) = 200 mjy, S(800) = 20 mjy and S(1250) 
= 2 mjy are shown. 



the rest-frame emission at wavelengths > 50/iin. Figure 8 shows the correla- 
tion between the luminosity in the starburst component, L s b, and the AGN 
dust torus component, L tor , for hyperluminous infrared galaxies, PG quasars 
(Rowan- Robinson 1995), and IRAS galaxies detected in all 4 bands (Rowan- 
Robinson and Crawford 1989) (this extends Fig 8 of Rowan- Robinson 1995). 
The range of the ratio between these quantities, with 0.1 < L sb /L tor < 10, 
is similar over a very wide range of infrared luminosity (5 orders of magni- 
tude), showing that the proposed separation into these two components for 
hyperluminous ir galaxies is not at all implausible. 

The radiative transfer models can be used to derive dust masses and hence, via 
an assumed gas-to-dust ratio, gas masses. For the M82-like starburst model 
used here the appropriate conversion is M^ust — 10~ 4 ' 6 L s fc, in solar units (Green 
and Rowan- Robinson 1996). Rowan- Robinson (2000a) has compared these es- 
timates with gas mass estimates derived from CO observations, where avail- 
able, and found good agreement (within a factor of 2). 
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8 10 12 14 16 

L sb 

Fig. 8. Bolometric luminosity in AGN dust torus component versus bolometric 
luminosity in starburst component: filled circles, hyperluminous ir galaxies (this 
paper); open triangles, PG quasars (Rowan- Robinson 1995); crosses, IRAS galaxies 
detected in 4 bands (Rowan-Robinson and Crawford 1989, galaxies with only upper 
limits on L tor omitted). 

The range of ratios of L s b/M gas for hyperluminous galaxies is consistent with 
that derived for ultraluminous starbursts. There is a tendency for the time- 
scale for gas- consumption, assuming a star formation rate given by eqn (1), to 
be shorter for the more luminous objects, in the range 10 7 — 10 8 yrs (alterna- 
tively this could indicate a higher value for the low-mass cutoff in the IMF). 
The cases where a strong limit can be set on M gas are also, generally, those 
where the seds do not support the presence of a starburst component. After 
correction for the effects of gravitational lensing, gas masses ranging up to 1-3 
x 1O 11 M are seen in most hyperluminous galaxies, comparable with the total 
stellar mass of an galaxy (10 1L2 (M/4L)/i5 2 ). Hughes et al (1997) argue 
that a star-forming galaxy can not be considered primeval unless it contains 
a total gas mass of 1O 12 M , but this seems to neglect the fact that 90 % of 
the mass of galaxies resides in the dark (probably non-baryonic) halo. 
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3 Conclusions and open issues 
cirrus 

* PAHs are a powerful diagnostic (uv surface brightness, metallicity) 

* need to consider illumination of interstellar dust by old (low-mass) and 
young (high-mass) stars 

* need to model the 3-D dust distribution (cf Silva et al 1998) 

* is the cold component far-out dust or giant grains ? 
starbursts 

* PAHS are a powerful diagnostic 

* need to model the evolution of a starburst (cf Efstathiou et al 2000) 

* need to consider escape of uv radiation durimg the high-r phase 

* need to consider phase when stars have dispersed the dust 
AGN dust tori 

* lack of 10 /iin feature a powerful constraint (need 10 \x spectroscopy) 

* need proper calculation of cloud model (ie 3-D radiative transfer) 

* how far out can dust torus extend ? A > 30/xm ? r > 350pc(L bo i / 10 13 ) ' 5 ? 

* is A em > 50/im emission from starburst ? Sub-mm interferometry could 
demonstrate the extended nature of the illuminating source. 

* there is a need for both an AGN dust torus and starburst components 
to understand most seds of hyperluminous ir galaxies. Measured gas masses 
support, in most cases, the starburst interpretation of rest-frame far-infrared 
and submm (A em > 50/zm)emission. 

Arp220-like starbursts 

* what fraction of ULIRGs are high-r starbursts ? This will be important for 
prospects for submm sepctroscopy with SCUBA and FIRST. 
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